A single chain analysis of doped quasi one dimensional 
spin 1 compounds: paramagnetic versus spin 1/2 doping 

3 



M. Fabrizio and R. Melin 



International School for Advanced Studies (SISSA-ISAS), 
Via Beirut 2-4, 34014 Trieste, Italy 
>• ■ and 

^ ! Istituto Nazionale per la Fisica della Materia (I.N.F.M.) 

\o 
o 
, 

Q\ , Abstract 

We present a numerical study of single chain models of doped spin 1 compounds. We 
use low energy effective one-dimensional models for both the cases of paramagnetic and 
spin-1/2 doping. In the case of paramagnetic doping, the effective model is equivalent 
to the bond disordered spin-1/2 chain model recently analyzed by means of real space 
renormalization group by Hyman and Yang. By means of exact diagonalizations in the 
XX limit, we confirm the stability of the Haldane phase for weak disorder. Above a critical 
amount of disorder, the effective model flows to the so called random singlet fixed point. In 
the case of spin-1/2 doping, we argue that the Haldane phase should be destabilized even 
for weak disorder. This picture is not in contradiction with existing experimental data. 
We also discuss the possible occurrence of (unobserved) antiferromagnetically ordered 
phases. 



1 Introduction 



Much effort has been devoted over the last decades to understand the behavior of quasi one- 
dimensional spin systems. These systems exhibit a rich variety of phases due to the enhanced 
quantum fluctuations in reduced dimensionality. One of the most famous and studied one dimen- 
sional property is the so called Peierls and spin-Peierls instability. Organic compounds exhibiting 
such an instability were synthesized already in the '70s [Q]. The more recent discovery of the inorganic 
CuGeOs spin-Peierls compound 0, ^] with a spin-Peierls temperature up to about 15 K renewed the 
interest to spin-Peierls systems, especially to the role of various possible doping mechanisms. In fact, 
it is possible to dope these inorganic compounds in a very controlled fashion with both magnetic 
and paramagnetic ions (see for instance [Q, || for Si doping, || [7|, |8| for Zn doping and || for Ni 
doping). It turns out that an antiferromagnetic (AF) phase is induced even for very small impurity 
concentrations, with a maximal Neel temperature of the order of 4 K at few percent doping. Neutron 
scattering experiments ||, ^ have unambiguously identified this phase as an antiferromagnet with 
a coexistence of gapless antiferromagnetic excitations at low energy and spin-Peierls excitations at 
higher energy H ||, and also with important disorder effects @. In a recent work, the present authors 
have shown that insight into the physics of these compounds can be gained by a single chain model 



[10, 11], which does exhibit an enhancement upon doping of low energy AF fluctuations, coexisting 
with higher energy spin-Peierls features. More specifically, by means of the renormalization proce- 



dure proposed in fl2|| , it was shown in [11] that this effective model belongs to the Random Singlet 
(RS) universality class of a random Heisenberg model [13]. Although within a single chain model it 
is difficult to handle rigorously the interchain couplings that stabilize the Neel phase, nevertheless 
this approach does emphasize the conjugated effects of low dimensional quantum fluctuations and 
disorder, which seems to be an important ingredient in the physics of these compounds. 

Another type of low dimensional quantum phase of spin systems is the Haldane gaped phase |14| 
of a spin-1 antiferromagnetic chain, which has been observed in quasi one dimensional compounds 
[|j~5|. Quite interestingly, a spin-1 chain and a dimerized spin- 1/2 chain belong to the same quantum 



phase, with a non vanishing string order parameter [16]. However, as we are going to argue, in spite 
of this analogy, the spin-Peierls and spin-1 chains behave very differently upon doping. Moreover, 
paramagnetic doping of spin-1 chains is very different from spin- 1/2 doping. In fact, there are some 
experimental evidences that the Haldane phase is robust against paramagnetic doping: it was found 



in [17] that the Haldane gap persists up to 20% doping in the NiC204DMIZ compound where Ni 
was substituted by Zn. On the other hand, it was shown in [18| that doping NENP compounds with 
less that 1% of Cu leads to a Curie behavior in the susceptibility, with a Curie constant 4.6 times 
larger than the one of the free Cu 2+ ions. These experiments thus show a very different behavior 
of quasi one-dimensional spin-1 compounds upon doping: weak paramagnetic doping preserves the 
Haldane gaped phase whereas weak spin-1/2 doping induces low lying energy states. We argue that 
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the behavior of quasi one dimensional spin-1 systems upon both spin-1/2 and paramagnetic doping 
can be understood on the same basis as the doped spin-Peierls compounds, namely, in the framework 
of a one dimensional effective model with a rigorous treatment of disorder effects. We also examine 
the issue of whether a strong doping antiferromagnetic phase could appear for paramagnetic doping. 
This article is organized as follows. In Section |2| we introduce an effective one dimensional model 



for paramagnetic doping. The zero temperature phases of this model were analyzed recently [19, 20 1. 
We check the nature of the zero temperature phases using exact diagonalizations of the effective 
spin-1/2 model at zero temperature in the XX limit. We then argue that, in close analogy to our 
approach for spin-Peierls compounds, the zero temperature Haldane phase of the effective model 
should be robust against switching interchain couplings whereas the RS phase should possibly turn 
into an antiferromagnetic phase. Section || is devoted to an analysis of an effective model for spin- 
1/2 doping. We show that in this case the Haldane phase is destabilized by an infinitesimal doping. 
In section |], we analyze the question whether an antiferromagnetic phase might appear for large 
paramagnetic doping. We carry out exact diagonalizations of small clusters, the interchain coupling 
being treated in mean field. These diagonalizations support our prediction, namely the possible 
emergence of a strong disordered antiferromagnetic phase. Finally, section ||| is devoted to some final 
remarks. 

2 Effective one dimensional model for paramagnetic doping 
2.1 The model 

Based on the Valence Bond (VB) description of a spin-1 chain Hyman and Yang [19] first 
argued that bond-disordered spin-1 chains can be described by the following low energy effective 
model: introducing bond disorder amounts as a first approximation to "break" the chain into clusters 
of consecutive strongly coupled spins. As a consequence, effective low energy spin-1/2 degrees of 



freedom appear at the edges of these segments [22 1 . Two neighboring 1/2-spins, one at the right and 



the other at the left edges of two consecutive segments are then assumed to be weakly AF coupled. 

This model should be also appropriate to describe the effects of doping by paramagnetic ions. 
In fact, if a magnetic ion is substituted by a non magnetic impurity (as in the case of Zn doping 
in NiC204DMIZ |l7|| ), one expects that two spin-1/2 edge excitations appear at the right and the 
left of the impurity, which get coupled by an AF exchange J2 weaker than the bulk exchange Jh- 
In addition, two spin-1/2 moments at the opposite edges of a segment are also coupled by virtual 
polarization of the spin-1 background. Assuming that these two edge moments Si and S2 are at 
distance /, the effective interaction mediated by the spin-1 background is [22] 

Wi, 2 (0 ~ J/K-lj'e-'/fcSi.Sa, (1) 

where we have discarded a \j\fl prefactor, irrelevant to the physics we want to discuss. 
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2.2 Exact diagonalizations in the XX limit 



Hyman and Yang |l9j, and also Monthus et al [20], carried out a RG analysis of this effective model. 
These authors showed that, for weak disorder, the system remains in the gaped Haldane phase. As 
the disorder increases above a critical value, the system exhibits a second order transition to the RS 
phase. The aim of this section is to check this prediction by means of exact diagonalizations in the 
XX limit. Due to the presence of ferromagnetic bonds, the XX limit is not expected to give a good 
description of the isotropic XXX model, contrary to what is found when all bonds are AF. However, 



like the XXX model analyzed in Ref. [19], also the XX version may have in principle two different 
phases for weak and strong doping. Therefore, if one is only interested in confirming the existence 
of this phase transition, the XX model should be sufficient. 

In the XX limit, this model can be maped onto non interacting spinless fermions via a Jordan- 
Wigner transformation p3|]: 

w = ^E J *( c H-i c *+ c i' c *+i)' (2) 

i 

where the label i denotes the sites carrying spin 1/2 moments and Jj is calculated according to the 
aforementioned rules. We restrict ourselves to an even number of spin 1/2 magnetic moments, so that 
the fermions are periodic and the Fermi sea is half-filled. Moreover, the impurities are distributed 
according to a Poisson law with an average impurity to impurity distance {I). 

Increasing the strength of disorder can be done either by decreasing J2 or by reducing the average 
impurity to impurity distance (I). We have found that in both cases increasing the strength of disorder 
leads to an increase of the correlation length and finally to power law correlations. We now present 
our calculations. 

In order to calculate the SqS^ correlations, we first diagonalize the tight binding Hamiltonian 

H (TB) = ^J i (\i + l)(i\ + \i)(i + l\). (3) 
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If \& Q = J2i N) denotes the eigenstates of (||), the average correlations are (S-Sj) = Aij — Bi/2 — 
Bj/2, with 

A i,j = E + E E*?*f*?*i (4) 

a,j3eFS a£FS /3&FS 

Bi = E W) 2 - (5) 

a£FS 

The correlations {SqS^) are shown on Fig. [l] in the strong disorder regime. This correlation function 
follows quite nicely a 1/R 2 decay, as in the random singlet phase [13]. 

We now keep a constant dilution of impurities but increase the nearest neighbor interaction J2 
across the impurities. This should drive the systems into the Haldane phase, with exponentially 
decaying correlation functions. This is clearly seen for the SqSr correlation plotted in Fig. || in a 
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10 100 1000 

Figure 1: Logarithm of the correlation between the z component of the spins in the strong disorder 
phase (paramagnetic doping), as a function of the spin-spin distance. The AF exchange between 
spins 1 is Jh = 1, the nearest neighbor AF exchange is J2 = 0.1, the correlation length in the 
Haldane phase is £0 = 20 lattice spacings, the average distance between two impurities is (/) = 50 
lattice spacings. For N = 100 (200) impurities, we performed the average over 32000 (2000) disorder 
realizations. The straight line is a fit to a l/R 2 behavior. 
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Figure 2: Logarithm of the correlation of the z component of the spins in the small disorder phase 
(paramagnetic doping). This is a semi-log plot. The AF exchange between spins 1 is Jh = 1; the 
nearest neighbor AF exchange is J2 = 0.5, 0.6, 0.7, 0.8; the correlation length in the Haldane phase 
is £0 = 20 lattice spacings; the average distance between two impurities is (I) = 50 lattice spacings. 
The correlations have been calculated for N = 100 impurities and averaged respectively over 50000, 
10000, 6000 and 1300 realizations of the disorder. The solid lines are a fit to the form exp (— R/£) /R 2 , 
with respectively £ = 23 ±0.5, 17.5 ±0.5, 14.5 ±0.5, 13 ±0.5. The saturation for correlations smaller 
than ~ 10~ 19 is due to round-off errors during the diagonalizations. 
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Figure 3: Logarithm of the correlation of the z component of the spins in the regime (/) < £o 
(paramagnetic doping). This is a log-log plot. The AF exchange between spins 1 is Jh = 1, the 
nearest neighbor AF exchange is J2 = 0.7 the correlation length in the Haldane phase is £0 = 20 
lattice spacings, the average distance between two impurities is (l) = 10, 15, 20 lattice spacings. The 
correlations have been calculated for N = 100 impurities and averaged respectively over 3400, 10000, 
25000 realizations of the disorder. The solid lines are a fit to the power law 1/R 2A in the case 
(I) = 10; to the form exp (— R/£)/R 2A with £ = 45 lattice spacing for (I) = 15; and £ = 30 lattice 
spacings for (I) = 20. 



semi- log plot for different values of J2. Notice on Fig.^ an increase of the scattering of the data as 
the system approaches the zero temperature critical point, in spite of an increase of the number of 
disorder realizations. This is a clear signature of critical fluctuations and thus consistent with the 



prediction that this zero temperature transition is a second order transition [19, 2C]. 

We now consider the case J2 = 0.7 on Fig. |2| and drive the system towards the infinite correlation 
length phase by decreasing the average impurity to impurity distance (/). The resulting correlations 
are plotted on Fig. ||. We clearly see that decreasing the average distance between impurities drives 
the system towards a phase with critical correlations. 

Although the effective model is valid only if (I) > £0, we have decided to investigate also the 
situation (I) < £0 just as a theoretical model. In this case we were not able to fit the correlation to 
the RS 1/R 2 form. Instead, we found in the case of Fig. || a power law decay different from the RS 
one. Indeed, we find 

1 



<«> ~ ( 6 ) 
This may be a finite size effect: we know that the renormalization group procedure is asymptotically 
correct since the disorder distribution becomes extremely large as the decimation is carried out. In 
the case (I) < £0, it is possible that, given the small sizes that we use, the distribution of bonds is 
not large enough, so that the correlation exponent has still not converged the one of the RS. Also, 
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we notice that the bond distribution is 

p (i J i)=2i|(S)" +%e <*-i J i) + ^ J -^ < 7 » 

In the case £o < (0; the smaller the exchange the larger the probability, which explains that in this 
regime the numerics works so well (see Figjl|). In the opposite regime (£q > (/)), the smaller the 
exchange the smaller the probability. In this case, we still expect the RS fixed point. However, it is 
not so obvious to exhibit it from a numerical point of view. 

3 Effective one dimensional model for spin 1/2 doping 

3.1 Effective one dimensional model 



The effect of spin 1/2 doping in quasi one dimension spin 1 NENP compounds was analyzed in [18] 



where susceptibility measurements were carried out, and also in [22] by means of ESR measurements. 
These last experiments showed unambiguously that spin 1/2 doping amounts to release two spin 1/2 
effective moments on the right and the left of the impurity, that are weakly AF coupled to the spin 
1/2 impurity. This picture is not in contradiction with the susceptibility measurements fl8fl , where 
the temperature variations of the susceptibility were found to have Curie behavior. The effective 
model thus consists of disordered units of three consecutive spin 1/2 spins. The spin 1/2 moments 
experience a weak AF interaction J' and two consecutive units are coupled via the ferromagnetic or 
AF interaction (g). 

It is immediately clear that the RG description of this effective model is very different from the 
one in The reason is that, in the weak disorder case, one first would replace all the 3 spin 

1/2 units by a spin 1/2. Therefore the effective model would be again a bond disordered S = 1/2 
Heisenberg AF, but the bonds being either ferromagnetic or AF. This model has been investigated 
in |Q] quite in detail and shown to be unstable against a small disorder. In particular, the ground 
state is expected to have power law correlations. We check in section |3.2| that our effective model 
indeed has power law correlations in its ground state. However, the physics is expected to be very 
different from the RS: the low temperature susceptibility has a Curie dependence whereas in the RS 
the low temperature susceptibility diverges slower than a Curie behavior (x(T) ~ l/(Tln 2 T)) []l3| . 
Moreover, the low temperature physics is dominated by large effective spins. Coming back to the 
physics of spin 1/2 doped quasi one dimensional spin 1 systems, we do not know whether, given this 
effective description, it would be possible to drive the system to a Neel phase by switching interchain 
couplings. 

3.2 Exact diagonalizations in the XX limit 



We repeat the numerical calculations of section |2.2] in the case of the effective model for spin 1/2 
doping. As in the case of paramagnetic doping, we only consider the spin anisotropic XX model. The 
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Figure 4: Logarithm of the correlations of the z component of the spins for the effective model of 
spin 1/2 doping. The data have been fited to a 1/R 2 dependence, we have chosen: Jh = 1, £o = 20, 
(I) = 50 and (a): J' = 0.8, 102 spin 1/2 moments (b): J' = 0.1, 102 spin 1/2 moments, (c): J' = 0.8, 
204 spin 1/2 moments, (d): J' = 0.1, 204 spin 1/2 moments. 

(SqSjj) correlations are plotted on Fig. |3] for weak and strong disorder and exhibit & 1/R 2 power law 
decay, like for the RS phase. This numerical calculation thus supports the previous argument that 
the effective one dimensional model suited for spin 1/2 doping would flow to a massless fixed point 
as soon as disorder is introduced. 



4 Exact diagonalizations of small clusters with paramagnetic dop- 
ing 

We now turn to exact diagonalizations of small clusters at finite temperatures, the interchain coupling 
being treated in mean field. Notice that this problem is quite difficult and we cannot go to quite large 
sizes because (i) one must calculate all the eigenvalues and eigenvectors since we are interested in the 
finite temperature behavior, (ii) one must calculate the staggered field in mean field, which implies 
solving for self consistent equations, (iii) one must average over disorder. This kind of calculation 



was already carried out in the context of disordered spin-Peierls systems [11], with results not in 
contradiction with experiments. We consider an eight site chain with two paramagnetic impurities, 
and thus six sites carrying a spin 1 in the squeezed chain. There are three inequivalent ways of putting 
the impurities on the chain as shown on Fig. |5[ We apply a finite staggered field and calculate the 
staggered magnetization response for different magnetic fields and temperatures. In practice, the 
temperature and the staggered field are varied between 0.01 and 2 with an increment of 0.01 and 
finally the self consistent staggered field h* is calculated as a function of temperature by imposing the 
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t— t — i .J* :: . (6) 

+ -4—^- t i > (6) 

• . i . > o) 

Figure 5: The three inequivalent disorder realizations, represented in the six sites squeezed chain. 
Arrows denote periodic boundary conditions. A dashed line denotes a weak bond. The sign of the 
staggered field is indicated by a ± sign. The degeneracy is indicated in brackets. 

condition h*{T) = J±m s (T, /i*(T)), where J± is the strength of the effective interchain coupling. The 
self consistent staggered field is plotted on Fig. || for two values of the effective interchain coupling 
J±. As in the case of the spin-Peierls system, we observe of Fig. ^ important fluctuations of the mean 
field Neel transition temperature. At this point, we conclude that a transition to an AF ordered state 
should be possible in the strong disorder regime. 

In order to have an idea of the variations of a macroscopic observable, we plotted on Fig.|?] the 
variations of the spin susceptibility as a function of temperature for two values of the interchain 
coupling. We observe that unlike the case of the spin-Peierls system , the emergence of AF as the 
temperature is lowered does not induce a maximum of the susceptibility, but a jump in the slope of 
the susceptibility versus temperature. 

Moreover, as a comparison between the two different doping mechanisms in a spin-1 chain, we 
note that eliminating the non magnetic sites going from the original chain to the squeezed chain 
has more drastic consequences for the paramagnetic doping than for the spin 1/2 doping. As far as 
paramagnetic doping is concerned, applying a staggered field in the original chain amounts to apply 
a staggered field in the same direction on both the spin 1/2 moments released by the impurity (on 
the right and the left of the impurity) and thus to systematically frustrate the AF coupling J2 across 
the impurity. Notice that in the case of spin 1/2 doping, this competition does not arise since the 
spin 1/2 impurity is coupled by AF bonds to the localized effective spin 1/2 moments. Hence, in the 
case of paramagnetic doping, this effect should play against the establishment of an AF phase. 

5 Conclusions 

We have thus carried out a detailed investigation of the effect of doping spin 1 chains by paramagnetic 
and spin 1/2 impurities. Our approach takes advantage of the quasi one dimensionality of these 
systems. We have shown how these systems can be maped onto bond disordered spin 1/2 chains, 
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Temperature 

Figure 6: Self consistent staggered field as a function of temperature calculated for Jh = 1, Ji = 0.5 
J± = -8 ((1),(2) and (3)) and J± = 0.4 ((4), (5), (6)). (1) and (4), (2) and (5), (3) and (6) correspond 
to the disorder realizations No 1, 2, 3 on Fig. [|. 
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Figure 7: Variations of the uniform susceptibility per spin as a function of temperature for Jh = 1, 
J 2 = 0.5 and (a): J± = 0.8 (b): J± = 0.4. 
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the zero temperature fixed points of which were recently analyzed (flJ^, ^] for the effective model 
appropriate to paramagnetic doping, and [24] for spin 1/2 doping). We have pointed out that 
paramagnetic and spin 1/2 doping are different physical situations since the effective one dimensional 
models have a very different physics. Weak paramagnetic doping preserves the Haldane phase [17|. 
Above a critical disorder strength, the effective one dimensional model flow to the RS phase. On 



the basis of our previous experience of spin-Peierls systems [10, 11], we argue that the RS phase, 
which was found to be the strong disorder zero temperature fixed point, would translate in these 
paramagnetically doped spin 1 chains into a possible AF phase above a critical disorder. Quasi one 



dimensional spin 1 compounds were already found to be robust against paramagnetic doping [17], 
but, to the best of the authors' knowledge, no AF ordered phase was found up to now. We have shown 
by means of exact diagonalizations of small clusters that a strong disorder AF phase is possible. In 
fact, the disorder favors the enhancement of magnetic fluctuations (as in the case of spin-Peierls 
compounds). Moreover, we have pointed out that in the paramagnetically doped spin 1 compounds, 
the interchain coupling always frustrates the AF coupling across the impurity. This mechanism does 
not exist in the spin 1/2 doped compounds, which in addition flow to a gapless phase for arbitrary 
weak disorder. 

Although the spin 1/2 doping seems at first sight more promising, it is not clear at this stage 
whether it should or not display an AF phase when interchain coupling is allowed, since the physics 
of the effective one dimensional model is dominated by clusters of spins coupled ferromagnetically, 
or, equivalently, by large spins. 

On the other hand, a very difficult problem in the possible observation of an AF phase in the 
paramagnetic doping case is that the microscopic parameters should be tuned in such way that AF 
occurs for not too large impurity concentrations, when the compounds become unstable. This means 
that the AF coupling across the impurity should be as small as possible. Therefore, the issue of 
the appearance of AF upon doping spin-1 quasi one dimensional compounds requires still further 
experimental investigations. 

Finally, even in the absence of AF ordering, neutron scattering experiments should reveal the 
presence of low-lying magnetic excitations along the chain directions, in both the spin-1 and spin-1/2 
doping situations. The zero temperature fluctuations should be critical for a small spin-1/2 doping. 
For spin-1 doping, the correlation length should remain finite below a critical amount of disorder, but 
the correlation length should increase with doping. To our knowledge, such measurements have not 
been carried out up to now in doped NENP or DMIZ compounds, even though neutron scattering 
experiments were performed in the doped quasi one dimensional transition oxyde Y2CaBaNii_. E Zn x .05 
[p5| . However, this compound is not properly described by the type of localized spin models that we 
have analyzed here. Nonetheless, the results in [25[] are interesting in view of our conclusions since 
the appearance of low lying AF fluctuations upon doping is reported. 
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